EFFECTIVE HAMILTONIANS FOR THIN DIRICHLET TUBES 
WITH VARYING CROSS-SECTION 



J. LAMPART, S. TEUFEL*, AND J. WACHSMUTH 



Abstract. We show how to translate recent results on effective Hamiltoni- 
ans for quantum systems constrained to a submanifold by a sharply peaked 
potential to quantum systems on thin Dirichlet tubes. While the structure 
of the problem and the form of the effective Hamiltonian stays the same, the 
difficulties in the proofs arc different. 

o 

> 
O 

| /' | The question whether a Schrddinger Hamiltonian, which localizes states close to 

a submanifold of the configuration space by large forces, may be replaced by an 
I effective operator on the submanifold is studied extensively and in various different 

settings in the literature (see e.g. [IJ|21II9fflG3ElEllHll§])- 
1 ^ 1 It is well-known that restricting the classical Hamiltonian system to the sub- 

manifold and then using Dirac's approach to quantizing constrained Hamiltonian 
I systems [TUj is too restricted. For there are lots of cases where the extrinsic cur- 

vature of the submanifold, which never shows up in Dirac's approach, plays a role. 
Therefore two other approaches have been investigated: 

£h • Soft constraints: A rapidly increasing potential is used to localize solutions 

close to the submanifold (see [TJ[2] and references therein). 
• Hard constraint: the localization is achieved via Dirichlet boundary con- 
ditions on a thin tube centered around the submanifold (see the reviews 

mm- 

If the potential or the tube's cross-section depend on the point on the submanifold, 
the constraint is called varying. First results for such constraints were given in 

Recently two of the authors have deduced effective Hamiltonians for the case of 
a varying soft constraint in arbitrary (co-)dimension (see [HI IS])- Here we explain 
how these results may also be obtained for the case of a varying hard constraint. 

1. The Setting 

Let (A, G) be a Riemannian manifold of dimension d + k and C a submanifold 
of dimension d without boundary, which is equipped with the induced metric g = 
G\tc- We assume that there is a non-self-intersecting tube Bs of radius S > 
around C. 

If C is compact, such a tube Bs always exists and is compact itself. For the sake 
of a simple presentation we will focus on the latter case in the following and only 
shortly comment on the necessary adjustments in the case of a non-compact C. 

Since there is a canonical diffeomorphism <£> from Bs into the normal bundle 
7r : NC — > C, we can scale any subset of Bs in the normal direction via 

D e :NC^NC, {q,N)^ (q,eN). 



Key words and phrases, thin tubes; effective Hamiltonians; constraints; spectral asymptotics. 

1 



EFFECTIVE HAMILTONIANS FOR THIN DIRICHLET TUBES WITH VARYING CROSS-SECTIOI2 

Let C B$ be an open subset with smooth boundary such that the cross-sections 

n{q) :=$(0)nAf,C 

are all diffeomorphic, compact, and connected. Then n : £1 — > C has the structure 
of a fiber bundle compatible with the one of NC. We assume that this bundle has 
smooth local trivializations. In the case of a non-compact C one has to postulate 
the existence of a set of local trivializations whose derivatives satisfy global bounds 
in a suitable manner. The e-thin tube O e is now defined via 

fl £ := Q^D^n. 

Our goal is to approximate the spectrum of and the unitary group generated by 

H e :=-e 2 A G on L 2 (ft e ,/x G ) 

with Dirichlet boundary conditions by using an effective Schrddinger operator 
on L 2 (C). Here A G is the Laplace-Beltrami operator associated with G and the 
factor e 2 has been put in for convenience because otherwise the spectrum of H e 
would diverge in the limit e —¥ 0. H e is obviously unitarily equivalent to the 
operator 

-£ 2 A$. G on L 2 ($(fi B ),fi*. ) 

with Dirichlet boundary conditions. We will identify the two operators in the 
following without making the diffeomorphism $ explicit anymore. 

2. Basic Ideas 

Consider the vector bundle £f := {(q, ip) | q G C, ip G C 00 (0(g))} over C, where 
the fibers f2(g) of the bundle fl are replaced with C°° (0(g)) and the bundle structure 
of fl is lifted by using the composition with the local trivializations of f2 as the new 
trivializations. Via the normal connection V on NC, which is induced by G, 
every vector t G T q C tangent to C can be lifted into the tangent spaces T q ^ n Q- of 
the corresponding fiber. The derivative of sections of £{ into the direction of the 
lift defines the so-called horizontal connection V h on £{ (see [5]). The associated 
Laplacian Ah coincides with the Laplace-Beltrami operator A g on C for functions 
that are constant on the fibers. 

As in [HI IH] the basic idea is that after a measure transformation and rescaling 
the normal coordinates n = N/e the Hamiltonian H e may be split as 

H £ = -e 2 A h -A n + 0{e). 

This suggests to define for each q G C the local fiber Hamiltonian 

H i {q):=-A n on L 2 (Vl(q),d\) 

with Dirichlet boundary conditions. Here dA is the Lebesgue measure induced from 
N q C ~ M fe . Since each fiber ft(q) is compact, the spectrum of Hf(q) is discrete for 
all q G C. At some fixed q we number the eigenvalues by J G No- Due to the 
smooth dependence of r2(g) on 1 this gives rise to continuous families of eigenvalues 
Ej(q), so-called energy bands. In general, these bands may cross. 

Definition 2.1. An energy band Ej is called admissible, if Ej(q) is simple for all 
q G C and the associated complex eigenspace bundle is trivializable, i.e., there is a 
global section ipj of normalized eigenf unctions. In addition, if C is non-compact, 
Ej has to satisfy a gap condition as in [Hj. 

As is well-known from the theory of elliptic operators, the lowest eigenvalue 
Eo(q) on the connected domain fl(q) is simple and po(q) can be chosen positive. So 
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the lowest energy band Eq is always admissible because the positivity of ipo ensures 
that it is a global section. For an admissible energy band Ej the subspace 

Vj := {^{x)yj{x,n) |V e L 2 (C,g)} C L 2 (Q,G) 

may be identified with L 2 (C,g) via the unitary operator 

Uj : Vj -> L 2 (C, g) , 1>(x)<pj(x, n) ^ ij>{x). 

V j is approximately invariant under H £ because the associated projector Pj satis- 
fies 

[Pj, H £ ] = [Pj, -e 2 A h ] + O(e) = O(e) 

in C(D(H £ ),'H). However, we are interested in the way the spectrum and the 
unitary group are affected by the geometry and the global structure of C. These 
effects are of order e 2 . Therefore we have to improve on the invariance of the 
subspaces. 



3. Results 

< oo. Via adiabatic perturbation theory it is possible to construct a 
projector Pj = Pj + ePj + e 2 P] and a unitary Uj-.V}^ L 2 (C) such that 

(1) [P £ j,H e ]x { -oo,E^]{H*) = 0{a% 



where X(-oo,B majl ] is the characteristic function of (—00, £ max ]. The construction of 
Pj is quite similar to the one in [8] . We comment on the differences below. Here we 
could in principle continue the construction to obtain a projector which is invariant 
up to errors of order e N for any N € N. 

Now we reformulate the main result from [5] for the case of thin Dirichlet tubes. 
Here we use the index formalism including the convention that one sums over 
repeated indices. Moreover, we use latin indices running from 1 to d for 

coordinates on C, greek indices a, /3, . . . running from d + 1 to d + k for the normal 
coordinates, and latin indices a, b, .. running from 1 to d + k for coordinates on Q. 

Theorem 3.1. Let Ej be an admissible energy band and £ max < 00. There are 
C < 00 and Eq > such that for all e < e$ there exist a Riemannian metric g e j 
on C, an orthogonal projection Pj, a unitary Uj : PjL 2 (Q, /iq) — > L (C,[i gj ) and 

H £ j := UjPjH £ PjUj* with domain UjD(H £ ), 

which satisfy the following: 

(a) Dynamics: Hj is self-adjoint on L 2 (C,fjL gj ) and 

f^-vr* - uf e~ iH ^ U}) Pj X (-oo,E max] (H e )\\ < Ce 3 \t\. 

(b) Spectrum: For all (E £ ) with limsup £ £ ,£ < E max one has 

(i) Hji/j e = E e ip e => \\(H £ ~ E £ )U £ j*ip £ \\ < Ce 3 \\U £ j*ip £ \\, 

(ii) H £ m £ = E £ W =>■ \\{H £ j - E £ )U £ jP £ j^ £ \\ < Ce 3 \\f £ \\. 

For ipi = X(-oo,.E max ] ( — e2 ^g + Ej)ipi the effective Hamiltonian Hj is given by 

(HHjiPik = J (gj lJ p!shP J ej A + W^Ax - e 2 ir 2 urR Hl (Ej)u £ 1 ip 1 

+ £ 2 V^ (Vgeom + Vbh + Kmb) Vl) d M S J + 0{s 3 ), 
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where 

+s^ al9 lm Wj Sm ( i pj\3n a n^ J ) n{qr 
P J S j = - e(^j|iV^j)o( g) - £ 2 K j ap(ipj\ln a n> :i id~ / ipj) n{q) 

+e 2 W^( yj | 2 (» a - ( ¥ > J |n« ¥ » J ))iV^j ) n(g)) 
R Ht {Ej) = {l-Pj)(H f -Ejy\l-Pj), 
U? = 2g^vf^~ J ed J +n a ^,W l Je 2 d 2 J , 

Vgeom — —^"Va + \H y — g a b + ^ aj + ^ y) i 

Vbh = ^'(V^KI-PjM^)^, 

Kmb = ^ 7 Q ^(5 7 ^j||n Q n ,3 9 5 v3,7) 0(g ), 

tm£/i W i/ie Weingarten mapping, n the mean curvature vector, 1Z and 1Z the Rie- 
mann tensors of C and A (see [S] for definitions of all the geometric objects). 

For a non-compact C additional bounds on the derivatives of tpj as in [5] are 
required. There a detailed discussion of the effective Hamiltonian is provided, too. 

As soon as ([T]) has been established, the proof of Theorem 3.1 goes exactly along 



the same lines as in [5]. The strategy to obtain ([!]) is also the same here, but the 



technical difficulties are different. The key facts that have to be derived are 

\c{D(Hr +1 ),D(H^)) 



(2) \\[-e 2 ^P.j]\\r (D(HT ^.D (H ^ = 0(e), 



(3) \\[-e z ^R Hi {E,j)]\\ c[D(HT+l)D{HT)) = O(e) 

for some I £ N and all m £ No- In addition, in [5] we had to make sure that the 
derivatives of (fj decay fast enough in the spatially infinite fibers and that their 
decay is not destroyed by application of energy cutoffs and resolvents, which is not 
necessary in the case considered here due to the boundary conditions. However, 
the boundary poses new problems in the proof of ^ k, On the one hand, the 
volume of the fibers is varying so that V h is only metric on sections which satisfy 
the Dirichlet condition. On the other hand, application of V h destroys the Dirichlct 
condition. Therefore one cannot only stick to the differential operators but has to 
make use of the spectral representation, too. 

Roughly speaking, ^ means to show that all the derivatives of tpj are uniformly 
bounded, in particular at the boundary. This can be done by locally mapping 
to the constant tube equipped with a suitable product metric and applying the 
procedures from [8]. Here the smoothness of the trivializations of fl enters. 

For (|]) one makes use of the fact that the fibers are compact so that the resolvent 
may be written as Rh { (Ej) — ~Y^j^jPi/{Ei — Ej). Then its derivatives may be 
controlled via some Weyl's law by choosing I large enough. 

4. Discussion of the Results 



Due to Theorem 3.1 the spectrum of H e is given, up to errors of order e 3 , by 
the spectra of Hj for Je No. With our approach it is possible to obtain not only 
the energies close to inf a(Hj) but also the excitations of order 1. In this energy 
regime the leading part —e 2 A g + Ej is a scmiclassical operator, whose dynamics 
explores distances of order 1 for times of order Therefore this is the relevant 
time scale, on which the global structure of an e-independent C is seen. Theorem 



3.1 allows to look at even much longer times. 
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The spectrum of — e 2 A g + Ej is quite well- understood. We discuss here the 
role of the corrections in Hj for constant Ej and an Ej with one non-degenerate 
minimum on a compact C. By using standard results from semiclassical analysis 
this discussion could be extended also to several degenarate minima. 

1) Ej is constant: In this case the level spacing of — e 2 A g + Ej is of order e 2 
close to Ej. So the low eigenvalues strongly depend both on the corrections in -p J € 
and on the effective potentials lg Com and Vj (see [9] for an example with global 
effects). Since the kinetic energy is small for eigenvalues close to Ej, the corrections 
to g and the off-band coupling Uf*RH f (Ej)Uf do not matter here. 

They only become relevant for energies of order 1 above Ej. Let < a < 2. 
According to Weyl's law, the level spacing at energies of order e a above Ej is of 
order e a 0-~ d / 2 )+ d m p or d = 1,2 this is always bigger than the approximation error of 
order e 3 . For d > 3 the approximation error is only smaller for a > (2d — 6) /(d — 2). 
Note that the minimal a is always strictly smaller than 2. To fully resolve the 
spectrum for energies of order 1 one would have to go to order d in the construction 
of the super-adiabatic projector and the effective Hamiltonian. However, even in 
cases where the effective Hamiltonian is not precise enough to resolve the small 
level spacing, Theorem |3.1| still yields good control over the dynamics of states in 
this energy regime on the relevant time scales. 

2) Ej has one non- degenerate minimum: Order 1 above supi?/ the level spacing 
of Hj is again given by Weyl's law, i.e., only of order e d . So we are in the same 
situation as for energies of order 1 above a constant Ej. Close to sup Ej no general 
statements can be made about the level spacing. Order 1 below sup Ej the spectrum 
of Hj is dominated by Ej resulting in a level spacing of order e. Thus the effective 
potentials, which are of order e 2 , may be ignored here. For energies of order 1 above 
inf Ej, however, the e-corrections to the kinetic energy become relevant because the 
eigenfunctions oscillate on a scale of order 

We sketch the emerging picture for the case of a closed curve C of length L: 




For energies above sup Ej the corresponding eigenfunctions are not localized but 
extended over the whole submanifold. Hence, global effects may occur here, too. 

5. Conclusions and Outlook 

We have derived an effective Hamiltonian for the problem of hard constraints 
in quantum mechanics that covers all the interesting energy scales. Although the 
technical difficulties in the soft and the hard constraint approach differ, the results 
have a very similar structure. This is due to the fact that the wave function con- 
centrates close to the submanifold. In the future, we will investigate the Laplacian 
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on thin Riemannian fiber bundles using our adiabatic techniques. To do so we split 
up the metric into a horizontal and a vertical part and scale only the latter by e: 

G = G\x ~\~ s G v . 

This is related to the so-called adiabatic limit in global analysis (see e.g. [H]). In 
this setting there is no concentration inside the fibers. Therefore an expansion of 
the metric must be replaced by an averaging procedure. As a consequence, the 
effective Hamiltonian will have a somewhat different structure. 
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